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Abstract 

Sphere packings in high dimensions have been the subject of recent theoretical interest. Em- 
ploying numerical and theoretical methods, we investigate the structural characteristics of random 
sequential addition (RSA) of congruent spheres in <i-dimensional Euclidean space M. d in the infinite- 
time or saturation limit for the first six space dimensions (1 < d < 6). Specifically, we determine 
the saturation density, pair correlation function, cumulative coordination number and the structure 
factor in each of these dimensions. We find that for 2 < d < 6, the saturation density (j) s scales 
with dimension as cj) s = a/2 d + c 2 d/2 d , where d = 0.202048 and c 2 = 0.973872. We also show 
analytically that the same density scaling persists in the high-dimensional limit, albeit with differ- 
ent coefficients. A byproduct of this high-dimensional analysis is a relatively sharp lower bound 
on the saturation density for any d given by S > (d + 2)(1 - S )/2 d+1 , where S G [0, 1] is the 
structure factor at k = (i.e., infinite-wavelength number variance) in the high-dimensional limit. 
We prove rigorously that a Palasti-like conjecture (the saturation density in R d is equal to that 
of the one-dimensional problem raised to the <ith power) cannot be true for RSA hyperspheres. 
We demonstrate that the structure factor S(k) must be analytic at k = and that RSA packings 
for 1 < d < 6 are nearly "hyperuniform." Consistent with the recent "decorrelation principle," we 
find that pair correlations markedly diminish as the space dimension increases up to six. We also 
obtain kissing (contact) number statistics for saturated RSA configurations on the surface of a d- 
dimensional sphere for dimensions 2 < d < 5 and compare to the maximal kissing numbers in these 
dimensions. We determine the structure factor exactly for the related "ghost" RSA packing in R d 
and show that its distance from "hyperuniformity" increases as the space dimension increases, ap- 
proaching a constant asymptotic value of 1/2. Our work has implications for the possible existence 
of disordered classical ground states for some continuous potentials in sufficiently high dimensions. 

PACS numbers: 05.20.-y, 61.20.-p 
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I. INTRODUCTION 



We call a collection of congruent spheres in <i-dimensional Euclidean space M. d a hard- 
sphere packing if no two spheres overlap. The study of the structure and macroscopic 
properties of hard spheres in physica 

a 



to at least the work of Boltzmann 



dimensions (d = 2 or 3) has a rich history, dating back 
"\. Hard-sphere packings have been used to mode 



variety of systems, including liquids |2j , amorphous and granular media [3] , and crystals "\ . 
There has been resurgent interest in hard-sphere packings in dimensions greater than three in 
both the physical and mathematical sciences. For example, it is known that the optimal way 
of sending digital signals over noisy channels corresponds to the densest sphere packing in a 

n 

high dimensional space |5j. These "error-correcting" codes underlie a variety of systems in 
digital communications and storage, including compact disks, cell phones and the Internet. 



Physicists have studied hard-sphere packings in high dimensions to gain insif 



it into ground 
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and glassy states of matter as well as phase behavior in lower dimensions [6, 

The determination of the densest packings in arbitrary dimension is a problem of long- 
standing interest in discrete geometry [5]. The packing density or simply density of a 
sphere packing is the fraction of space M. d covered by the spheres, i.e., 

<P = pv 1 (R), (1) 

where p is the number density, 

n d/2 

V JR) = — —R d (2) 

u ; r(l + d/2) v ; 

is the volume of a ci-dimensional sphere of radius R, and T(x) is the gamma function. We 
call 

0max = SUp </>(P) (3) 



the maximal density, where the supremum is taken over all packings in M. d . The sphere 
packing problem seeks to answer the following question: Among all packings of congruent 
spheres, what is the maximal packing density </> max , i-e., largest fraction of WL d covered by the 
spheres, and what are the corresponding arrangements of the spheres (5, n|? For d — 1,2 
and 3, the optimal solutions are known. For 3 < d < 10, the densest known packings of 
congruent spheres are Bravais lattice packings [5|], but in sufficiently large dimensions the 
optimal packings are likely to be non-Bravais lattice packings. Upper and lower bounds on 



the maximal density max exist in all dimensions [5j. For example, Minkowski 11] proved 
that the maximal density ^ H among all Bravais lattice packings for d > 2 satisfies the 
lower bound 



^ >M ( 4 ) 

max — 2^~1 ' 



where ((d) = Y^=i k~ d * s the Riemann zeta function. It is seen that for large values of d, 
the asymptotic behavior of the nonconstructive Minkowski lower bound is controlled by 2~ d . 
Note that the density of a saturated packing of congruent spheres in M. d for all d satisfies 

which has the same dominant exponential term as (@J). A saturated packing of congruent 
spheres of unit diameter and density <fi in 3t d has the property that each point in space 
lies within a unit distance from the center of some sphere. In the large-dimensional limit, 
Kabatiansky and Levenshtein [3] showed that the maximal density is bounded from above 
according to the asymptotic upper bound 

1 



'max 



— 20-5990d' 



The present pap er is motivated by some recent work on disordered sphere-packings in high 
dimensions QQ. In Ref. Q, we introduced a generalization of the well-known random 
sequential addition (RSA) process for hard spheres in (^-dimensional Euclidean space M. d . 
This model can be viewed as a special "thinning" of a Poisson point process such that the 
subset of points at the end of the thinning process corresponds to a sphere packing. One 
obvious rule is to retain a test sphere at time t only if it does not overlap a sphere that 
was successfully added to the packing at an earlier time. This criterion defines the standard 
RSA process in M. d ^, , which generates a homogeneous and isotropic sphere packing in 
M. d with a time-dependent density (pit). In the limit t — > oo, the RSA process corresponds 
to a saturated packing with a maximal or saturation density 0(oo) = lim^oo <p(t) [16]. 
In one dimension, the RSA process is commonly known as the "car parking problem", 



which Rehyi showed has a saturation density 0(oo) = 0.747598. . . 17J. For 2 < d < oo, 
an exact determination of (t>(oo) is not possible, but estimates for it have been obtained 

fin 

via computer experiments in two dimensions (circular disks) |18l . I20j and three dimensions 
(spheres) j^, 22 1. However, estimates of the saturation density 0(oo) in higher dimensions 
have heretofore not been obtained. 
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Another thinning criterion retains a test sphere centered at position r at time t if no 
other test sphere is within a unit radial distance from r for the time interval nt prior to t, 
where k is a positive constant in the interval [0, 1]. This packing is a subset of the RSA 
packing, which we call the generalized RSA process. Note that when k = 0, the standard 
RSA process is recovered, and when k = 1, we obtain the "ghost" RSA process Q], which is 
amenable to exact analysis. In particular, we showed that the n-particle correlation function 
<7n(ri, r 2 , . . . , r n ) for the ghost RSA packing can be obtained analytically for any n, all 
allowable densities and in any dimension. This represents the first exactly solvable disordered 
sphere-packing model in arbitrary dimension. For statistically homogeneous packings in M. d , 
these correlation functions are defined so that p n g n (ri, r 2 , . . . , r n ) is proportional to the 
probability density for simultaneously finding n particles at locations r 1; r 2 , . . . ,r n within 
the system, where p is the number density. Thus, in a packing without long-range order, 
each g n approaches unity when all particle positions become widely separated within M. d , 
indicating no spatial correlations. The fact that the maximal density 4>(oo) = l/2 d of the 
ghost RSA packing implies that there may be disordered sphere packings in sufficiently high 
d whose density exceeds Minkowski's lower bound ©. 

n 

Indeed, in Ref. [14], a conjectural lower bound on the density of disordered sphere pack- 
ings Q| was employed to provide the putative exponential improvement on Minkowski's 
100-year-old bound. The asymptotic behavior of the conjectural lower bound is controlled 
by 2~(°' 77865 ---) d . These results suggest that the densest packings in sufficiently high dimen- 
sions may be disordered rather than periodic, implying the existence of disordered classical 
ground states for some continuous potentials. In addition, a decorrelation principle for disor- 
dered packings was identified ^| , which states that unconstrained correlations in disordered 
sphere packings vanish asymptotically in high dimensions and that the g n for any n > 3 
can be inferred entirely from a knowledge of the number density p and the pair correlation 
function g^iy)- At first glace, one might be tempted to conclude that the decorrelation 
principle is an expected "mean-field" behavior, which is not the case. For example, it is 
well known that in some spin systems correlations vanish in the limit d — > oo and the sys- 
tem approaches the mean-field behavior. While this notion is meaningful for spin systems 
with attractive interactions, it is not for hard-core systems. The latter is characterized by 
a total potential energy that is either zero or infinite, and thus cannot be characterized by 
a mean field. Furthermore, mean-field theories are limited to equilibrium considerations, 
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and thus do not distinguish between "constrained" and "unconstrained" correlations that 
arise in non-equilibrium packings of which there an infinite number of distinct ensembles. 
The decorrelation principle is a statement about any disordered packing, equilibrium or 
not. For example, contact delta functions (constrained correlations) are an important at- 
tribute of non-equilibrium jammed disordered packings and have no analog in equilibrium 
lattice models of any dimension. Finally, the decorrelation principle arises from the fact 
that the Kabatiansky-Levenshtein asymptotic upper bound on the maximal packing density 
(jHJ) implies that must go to zero at least as fast as 2~ 5990d for large d and therefore, 
unconstrained spatial correlations between spheres are expected to vanish, i.e., statistical 
independence is established ^j. There is no counterpart of the Kabatiansky-Levenshtein 
bound in mean-field theories. 

Motivated by these recent results, we study the structural properties of standard RSA 
packings (k = for the generalized RSA process) in the infinite-time or saturation limit, 
generated via computer simulations, for the first six Euclidean space dimensions (1 < d < 6) . 

nrrn 

The algorithm is checked by reproducing some known results for d = 1, 2 and 3 [13, UM, 20, 
l2lL I22B . Although we know that the saturation density (f>(oo) is bounded from below by 2~ d 
|l3l ll4||. the manner in which 0(oo) scales with density is not known for d > 3. One objective 
of this paper is to answer this question. Another aim is determine the corresponding pair 
correlation functions and structure factors in order to ascertain whether decorrelations can 
be observed as the space dimension increases up to six. A byproduct of our high-dimensional 
analysis is a relatively sharp lower bound on the saturation density of RSA packings for any 
d. Although a Palasti-like conjecture (the saturation density in M. d is equal to that of the 
one-dimensional problem raised to the dth power) is exact for ghost RSA packings, we 
rigorously prove that this conjecture cannot be true for standard RSA packings. 

In Appendix A, we obtain kissing (contact) number statistics for saturated RSA configu- 
rations on the surface of a (i-dimensional sphere for dimensions 2 < d < 5 and compare to the 
maximal kissing numbers in these dimensions. In Appendix B, we determine the structure 
actor exactly for ghost RSA packings and show that its distance from "hyperuniformity" 
23| increases as the space dimension increases, approaching a constant asymptotic value. 
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II. SOME KNOWN ASYMPTOTIC RESULTS FOR RSA PACKINGS 



Here we collect some known asymptotic results for the standard RSA process for d- 
dimensional hard spheres. Henceforth, we call <ft s = <f){oo) the saturation (infinite-time) 
limit of the density. 

In his numerical study of RSA hard disks, Feder postulated that the asymptotic 
coverage in the long-time limit for <i-dimensional hard spheres follows the algebraic behavior 

<P S - 0(r) ~ r-^ d , (7) 

where r represents a dimensionless time. Theoretical arg uments supporting Feder's law (JJJ) 



have been put forth by Pomeau |24| and Swendsen |25|. Not surprisingly, the saturation 
limit is approached more slowly as the space dimension increases. 

Moreover, similar arguments lead to the conclusion that the pair correlation function 
<72(f) at the saturation limit possesses a logarithmic singularity as the dimensionless radial 
distance r for spheres of diameter D approaches the contact value, independent of dimension 



24, 



25j, i.e., 

<?2(r) ~ ln(r — D), r — > D and = <p s . (8) 



Boyer et al. 26] also showed that the pair correlation function for d = 1 has super- exponential 
decay. Specifically, they found that at any finite time r or density 0, 

If 2 X'~ D 
9 ^ r ^V{r)\ hx{r-D) ) ' r ^ 00 and < ^ < ^- ( 9 ) 

Thus, $2 is a short-ranged function at any density. This super-exponential decay of the pair 
correlation function persists in higher dimensions as well. As we will discuss in Section IV, 
this rapid decay of (r) has implications for the analytic properties of the structure factor 
S(k). 

III. NUMERICAL PROCEDURES 

In what follows, we describe an efficient procedure to generate RSA packings in the 
saturation limit as well as the methods used to compute structural information, such as the 
density, pair correlation function, structure factor and cumulative coordination number. 
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A. Generation of RSA Packings in M. d in the Saturation Limit 



We present a computationally fast method to generate RSA configurations of hard spheres in 
the saturation limit in the thermodynamic limit. Periodic boundary conditions are applied 
to a hypercubic fundamental cell of side length L and volume L d . Spheres of diameter D are 
placed randomly and sequentially inside the fundamental cell, which is periodically replicated 
to fill all of d- dimensional Euclidean space M. d , until the saturation limit is achieved. 

In order to speed up the computation, we attempt to add a particle only in the available 
space rather than wasting computational time in attempting to add particles anywhere in the 
fundamental cell 27, 3|. This requires keeping track of the time- dependent available space, 
which is the space exterior to the union of the exclusion spheres of radius D centered at each 
successfully added sphere at any particular time. This is done by tessellating the hypercubic 
fundamental cell into smaller, disjoint hypercubic "voxels," which have side length between 
0.025-D and 0.1.D, depending upon the dimension. 

At the start of the simulation, all voxels are declared accessible to particle placement. 
There are two stages involved to determine the time-dependent available space. A coarse 
estimation of the available space is used in the first stage, which is refined in the second 
stage. In the first stage, a particle is successfully added to the simulation box, provided 
that it does not overlap any existing particle. To avoid checking for nonoverlaps with every 
successfully added particle to the simulation box, we employ a "neighbor list" |3[, which 
amounts to checking within a local neighborhood of the attempted particle placement. For 
each successfully added particle, all voxels located within the largest possible inscribed 
hypercube centered at the exclusion sphere of radius D that are fully occupied by the 
particle are declared to be part of the unavailable space. The voxels outside this hypercube 
but within the exclusion sphere may be partially filled. In the initial stages, such partially 
filled voxels are declared to be part of the available space. We call these accessible voxels. If 
there have been at least one million unsuccessful placement attempts since the last accepted 
particle placement, we move to the second stage to refine our determination of the available 
space. In particular, we determine whether each remaining accessible voxel from the first 
stage can accommodate a particle center by a random search of each accessible voxel. After 
about 1000 random placement attempts, a particle is either added to a particular voxel or 
this voxel is declared to be part of the unavailable space. This search is carried out for 
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all other accessible voxels. The simulation terminates when all voxels are unavailable for 
particle addition in the second stage. 

This two-stage procedure enables us to generate RSA packings that are saturated or 
nearly saturated. Generating truly saturated RSA packings becomes increasingly difficult 
as the space dimension increases, as the asymptotic relation Q indicates. 



B. Calculation of the Saturation Density 

At any instant of time r, the number N(t) of added particles for a particular configuration 
is known and the density </>(t) is computed from relation (0) with p = N{r)/L d . We call 
4>sto P = §{j max ) the "stopping" density, i.e., the density at the time r max when the simulation 
is terminated. The system size L/D is sufficiently large so as produce a histogram for <p ma x 
or (p s that is Gaussian distributed. Although we do not present the full distribution of 
densities here, we do report the associated standard errors. In order to estimate the true 
saturation density <f) s , the volume of the available space V a {r) as a function of dimensionless 
time t is recorded in the very late stages, namely, for the last ten particles added. The 
saturation density <p s is estimated from this late-stage data by plotting 0(r) versus r _1 / d 
[cf. (J7J)] and extrapolating to the infinite-time limit. However, to perform the extrapolation 
properly the time increment between each particle addition cannot be taken to be uniform 
but instead must increase with increasing time in order to account for the fact that we only 
attempt to add particles in the available space. In the very late stages, this time increment 
At is given by 

AT ^W (10) 

The stopping density (p s to P always bounds the saturation density (fi s from below, but as we 
will soon see, stop is very nearly equal to the saturation density (f> s . 



C. Calculation of the Pair Correlation Function 



We obtain the pair correlation function g^O") at the nearly-saturated stopping density 
4>stop for a specific configuration by generating a histogram of the average number of particle 
centers n(r) contained in a concentric shell of finite thickness Ar at radial distance r from 
an arbitrary reference particle center (^J . The radial distance r is defined as halfway between 
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the inner radius (r — Ar/2) and the outer radius (r + Ar/2) of each shell. The shell thickness 
is termed the bin width. Let rifc(r) represent the accumulated pairs of particles for the entire 
system placed in bin k associated with a radial distance r. By definition, n^(r) must be an 
even integer. Then 

n (r) = —jj- , (11) 

where N is the number of particles in the fundamental cell. In general, the pair correlation 
(or radial distribution) function is defined as 

92(r) = , (12) 

pv s heii{r) 

where v s h e ii is the volume of the <i-dimensional shell, given by 

(r + Ar/2) d - (r - Ar/2) dl 



V shell = «i(r) 



yd 



(13) 



p is the number density N/L d , and V\(r) is the volume of a rf-dimensional sphere of radius 
r as shown earlier. 

We compute ensemble-averaged pair correlation functions by binning up to a maximum 
distance of r max for each realization of the ensemble and then averaging over all ensemble 
members. Away from contact, we employ a bin width of Ar = 0.05.D. Near contact, we use 
a finer bin width of Ar = 0.005.D in order to accurately capture the logarithmic divergence 
of #2( r ) a s the contact value is approached. 



D. Calculation of the Cumulative Coordination Number 



Another quantity of interest is the cumulative coordination number Z(r), which gives 
the average number of sphere centers within a distance r from a given sphere center. It is 
related to the pair correlation function g 2 (r/D) as follows: 

/r/D 
Si(x)g2(x)dx 

/r/D 
x d - l g 2 (x)dx, (14) 

where s\(r) = d7c d ^ 2 r d ^ 1 /T(l + d/2) is the d- dimensional surface area of a sphere of radius 
r 0. 
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E. Calculation of the Structure Factor 



Finally, we also compute the structure factor S(k), which provides a measure of the 
density fluctuations at a particular wave vector k and is defined by the relation 

S(k) = l + ph(k), (15) 

where h(k) is the Fourier transform of the total correlation function h(r) = g 2 (r) — 1. When 
the total correlation in M. d depends on the wavenumber k = |k|, then the structure factor 
S(k) in M. d for any space dimension d is given by 



d 



OO 



S(k) = l + p (27^ / r^h{r) J ^;^} dr, (16) 



lo ' ' {kr) {d/2) - 1 

where J u (x) is the Bessel function of order v. 

The expression ()16|) provides a means for computing the structure factor by Fourier 
transforming the real-space total correlation function in M. d . If one is interested in the large- 
wavelength (small k) behavior, however, the large r behavior of h(r) must be known with 
high precision. Even for relatively large simulation cells, it is difficult to access this large-r 
asymptotic behavior. In such instances, it is better to compute the structure factor directly 
from the collective density variables, i.e., 

*<k) = ^, (17) 

where 

N 

p(k) = ^exp(zk-r J ) (18) 
i=i 

are the collective density variables, angular brackets denote an ensemble average, and k 
are the wave vectors appropriate for the periodic cell of volume V. For the hypercubic 
fundamental cell of side length L considered here, the d- dimensional wave vectors are given 
by 

/2tt 2tt 2tt \ 

k = I — ni, — n 2 , ■ ■ ■ , — n d 1 , (19) 

where (i — 1, 2, . . . , d) are the integers. Thus, the smallest positive wave vector that one 
can measure has magnitude 2n/L. For small to intermediate values of k, we will employ the 
direct method, while for intermediate to large values of k, we will use both the direct and 
indirect method [i.e., we calculate S (k) using (flfijl ]. 
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TABLE I: The computed saturation density <p s and associated standard error for the first six space 
dimensions. Included in the table is the stopping density 4> s top, relative system volume L d /v\{l/2), 
and the total number of configurations n con f. 



Dimension, d 


4^ stop 




L d /vi (1/2) 


' v conj 


1 


0.74750 


0.74750 ± 0.000078 


6688.068486 


1000 


2 


0.54689 


0.54700 ± 0.000063 


9195.402299 


1000 


3 


0.38118 


0.38278 ± 0.000046 


13333.333333 


1000 


4 


0.25318 


0.25454 ± 0.000091 


21390.374000 


635 


5 


0.16046 


0.16102 ± 0.000036 


66666.666667 


150 


6 


0.09371 


0.09394 ± 0.000048 


193509.198363 


75 



IV. RESULTS AND DISCUSSION 



A. Saturation Density 

The saturation density <p s for each of the first six space dimensions was determined by 
considering 100-1,000 realizations and several different system sizes, as discussed in Section 
III.B. Table U] summarizes our results for the saturation density for the largest systems 
and the associated standard error. Included in the table is the stopping density 4> s top, 
relative system volume L d /vi(l/2) for the largest system, where Vi(l/2) is the volume of 



a hypersphere, and the total number of configuration s n rqn f. The results for d = 1, 2 and 
3 agree well with known results for these dimensions 



iSJWl TI i 



21 



22 



. We see that the 



stopping density stop is very nearly equal to the saturation density (p s for all dimensions, 
except for d — 1 where these two quantities are identical. For d — 1, no extrapolation was 
required since we can ensure that the packings were truly saturated in this instance. 

It is of interest to determine how the saturation density <ft s scales with dimension. We 
already noted that the infinite-time density of the ghost RSA packing (equal to 2 d ) provides 
a lower bound on saturation density of the the standard RSA packing. Therefore, it is 
natural to consider the ratio of the saturation density to the infinite-time density of the 
ghost RSA packing, i.e., 2 d cf) s . When this ratio is plotted versus dimension for 2 < d < 6, 
it is clear that the resulting function, to an excellent approximation, is linear in d, implying 
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the scaling form 



£1 c 2 d 
2 d 2 d ' 



(20) 



where c\ = 0.202048 and c 2 = 0.973872. Indeed, the linear fit of 2 d (f) s , shown in Fig. d is 
essentially perfect (the correlation coefficient is 0.9993). If the data shown in Fig. ^is instead 
fitted up to d = 5, the predicted density for d = 6 from this fit function is within 0.4% of 
the corresponding datum. This indicates that the scaling form for relatively low dimensions 
is accurate. In the following subsection, we provide an analytical argument supporting the 
same scaling form in the high- dimensional limit. It is noteworthy that the best rigorous 
lower bound on the maximal density [^flj . derived by considering lattice packings, has the 
same form as (j2~Uj) . 

7 




FIG. 1: (Color online) Fit of data for the product 2 d (f) s to the linear form 1)20(1 for 2 < d < 6. The 
correlation coefficient is 0.9993, and a = 0.202048 and c 2 = 0.973872. 

An interesting conjecture due to Palasti [3^ claims that the saturation density for RSA 
packings of congruent, oriented d-dimensional cubes equals the saturation density S = 
0.747598 . . . 0] of the one-dimensional problem raised to the c?th power. It took over thirty 
years to show, through a precise Monte Carlo simulation in two dimensions 27j], that the 
Palasti conjecture could not be rigorously true. For RSA packings of congruent, oriented 
squares, the saturation density was determined to be 0.562009 ± 0.000004, which is close 
but not equal to (0.747598 . . .) 2 = 0.5589 .... It is noteworthy that the saturation density 
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for oriented squares is also close to that of circular disks (see Table |TJ). These two systems 
are distinguished from one another in that the available space for particle addition in the 



late stages for the former are relatively large rectangles, 27( while for disks they are small 
"triangular" -shaped regions j^ol ]. 

Is a Palasti-like conjecture (involving raising the one- dimensional density result to the dth 
power) ever valid for disks? We make the simple observation here that the Palasti conjecture 
is exact for the ghost RSA packing in the infinite-time limit because <f>(oo) = 2~ d in any 
dimension. Moreover, it is trivial for us to rigorously prove that a Palasti-like conjecture 
cannot be true for the standard RSA packing of spheres in M d . This conjecture would state 
that 

^ s = 2(0.419665. ..)d (21) 

is the saturation density for such a packing for all d. However, this violates the asymptotic 
Kabatiansky-Levenshtein upper bound (jHJ) for the maximal density of a sphere packing in 
M. d . Therefore, a Palasti-like conjecture cannot be true for standard RSA packing of spheres 
M. d . This was known from numerical experiments, but a proof was never presented until 
now. 



B. Pair Correlation Function 

Figures El and |3] show the ensemble averaged pair correlation functions for the first six 
space dimensions very near their respective saturation densities. They are computed from 
the same configurations used to calculate the saturation densities, as described in Section 
IIV Al To our knowledge, our results for $2 very near the saturation densities have not been 
presented before for d > 3. The inset in each figure shows the near-contact behavior, which 
is consistent with the expected logarithmic divergence at contact and fitted to the form: 

g 2 (x) =a ln(;r-l)+ai, 1 < x < 1.135, (22) 

where x = r/D. Table ITD summarizes the values of the fit parameters a Q and a\ for each 
dimension. Of course, the logarithmic term overwhelms the constant coefficient ci\ as x — i . 
Our result for the logarithmic coefficient ao for d — 1 agrees well with the exact result 
a Q = —1.128 . . . 20]. There are no exact results for an for d > 2, but it has been previously 
evaluated numerically for d = 2 by Hinrichsen et al |20(, who obtained the value a = —1.18, 
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TABLE II: Fits of the pair correlation function g<i (r) at the very nearly saturation density <fi = (p s top 
to the form ao ln(:c — 1) + a% for 1 < x < 1.135. 



Dimension, d 






1 


-1.11955 


-0.117475 


2 


-1.29175 


-0.883021 


3 


-1.16546 


-0.808843 


4 


-1.00743 


-0.58044 


5 


-0.714731 


0.020810 


6 


-0.412100 


0.661348 



which is somewhat smaller in magnitude than the value reported in Table ITT1 These authors 
were only able to fit their data in the near-contact region over about 1.5 decades on semi- 
logarithmic plot due to insufficient statistics. Indeed, we have employed substantially more 
configurations than they did and were able to fit our data over about 4.6 decades on semi- 
logarithmic plot. The results reported in Table ILT1 for d > 3 have not been presented before. 
Feder et al. also gave an expression for the dominant logarithmic term for any d in terms 
of a certain Voronoi statistic and the "hole-size" distribution function at contact; but since 
neither of these quantities are known analytically, it is not a practically useful relationship. 

Figure 0] plots all of the pair correlation functions on the same scale. We see that the 
"decorrelation principle," which states that unconstrained spatial correlations diminish as 
the dimension increases and vanish entirely in the limit d — > oo jl^l . Il^ , is already markedly 
apparent in these relatively low dimensions. Correlations away from contact are clearly 
decreasing as d increases from d = 1. The near-contact behavior also is consistent with the 
decorrelation principle for d > 2. Although the logarithmic coefficient ao increases in going 
from d = 1 to d = 2, it decreases for all d > 2 [cf. Table ILT] . The decorrelation principle 
dictates that ao tends to zero as d tends to infinity. Similarly, the constant coefficient a% 
increases for d > 2 and for d = 6 is equal to 0.661348 [cf. Table ILT] . Indeed, the decorrelation 
principle requires that a\ tends to unity, indicating the absence of spatial correlations, as d 
tends to infinity. 
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FIG. 2: (Color online) The pair correlation functions for RSA packings for d = 1 (top panel), 
d = 2 (middle panel) and d = 3 (bottom panel) very near their respective saturation densities: 
4> = <frstop = 0.74750, cp = 4> s top = 0.54689 and cp = 4> s top = 0.38118, respectively. The insets show 
semi-logarithmic plots of the divergence in g2 (r) near contact. The straight line is a linear fit of 
the data. 



C. Cumulative Coordination Number 



The cumulative coordination number Z(r) at <fi = <p s top is easily obtained from the previ- 
ous results for (?2( r ) by performing the integration indicated in (fT^j) . All of these results for 
1 < d < 6, to our knowledge, have not been presented before. For D < r < D(l + e), where 



16 




FIG. 3: (Color online) The pair correlation functions for RSA packings for d = 4 (top panel), 
d = 5 (middle panel) and d = 6 (bottom panel) very near their respective saturation densities: 
4> = <frstop = 0.25318, cp = 4> s top = 0.16046 and cp = 4> s top = 0.09371, respectively. The insets show 
semilogarithmic plots of the divergence in Q2 (r) near contact. The straight line is a linear fit of the 
data. 



e — > 0, we can use the asymptotic form to yield the corresponding expression for Z(l + e) 
in any dimension d as 



Z(l + e) = T 



a {ln(e) — 1} + a x 



e + o e 



(23) 
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FIG. 4: (Color online) The pair correlation functions for RSA packings for the first six dimensions 
very near their respective saturation densities. Correlations clearly decrease as the space dimension 
increases. Note that the first intercept of g<i (r) with unity decreases with increasing dimension. 

where o(e) signifies terms of higher order in e and distance is measured in units of the 
hard-sphere diameter. Thus, Z(l + e) — > in the limit e —>■ 0, i.e., the average contact 
number is zero for RSA packings. This feature makes RSA packings distinctly different 
from maximally random jammed (MRJ) packings 0, [3| , which have an average contact 
number equal to 2d j3, 13 ■ This is one reason, among others, why the former packing 
has a substantially smaller density than the latter. The MRJ densities, as determined from 
computer simulations 



31 



32 



33, 



34j, are given 0.64, 0.46, 0.31 and 0.20 for d = 3,4,5 and 



6, respectively, which should be compared to the RSA saturation densities given in Table |TJ 
We also note that the appearance of the product ln(e)e in ()23|) means that the cumulative 
coordination number will be concave near contact and possess a positive infinite slope at 
contact. 

For values of r away from the near-contact behavior, we can deduce the following ex- 
plicit approximation for the <i-dimensional RSA cumulative coordination number Z(x) as a 
function of the dimensionless distance x = r/D at saturation: 



Z(x) = (c\d + c 2 d 



a {m(e) — 1} + ai 



x 



d 



x > 1 + e, 



(24) 



where e is a small positive number (which can be taken to be 0.135 in practice) describing the 
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range of the near-contact behavior, and the constants ci,C2,o,Q,ai are given in the caption 
of Fig. ^ and Table |Hj This approximation is obtained using the saturation density scaling 
(I2(J|) . the near-contact relation (J23|) . and definition (|14j) employing the approximation that 
92{x) — 1 ; which of course becomes exact as x and/or d becomes large. In light of the 
superexponential decay of #2 m any dimension and the decorrelation principle, x or d does 
not have to be large for the approximation ()24|) to be accurate. 

Figure El shows the cumulative coordination number Z(r) for the first six space dimensions 
at their respective saturation densities. These results are obtained by numerically integrating 
(jl4)l using the trapezoidal rule and our corresponding numerical data for g^O")- The insets 
of these figures clearly show the concavity of Z(r) near contact, as predicted by (|23|) . Away 
from contact, formula (J24j) provides a good approximation to the numerically determined 
values of Z[f) and is especially accurate for d > 3. 

In Appendix A, we determine kissing (contact) number statistics for saturated RSA 
configurations on the surface of a <i-dimensional sphere for dimensions 2 < d < 5 and 
compare to the maximal kissing numbers in these dimensions. It is of interest to determine 
the value of r at which the cumulative coordination number Z(r) for a saturated RSA packing 
in lBL d matches the average RSA kissing number (Z) on the surface of a hypersphere in the 
same dimension. Using Table IPTl we find that Z(r) = (Z) for r/D = 1.4233, 1.43042, 1.36044 
and 1.3202 for d = 2,3,4 and 5, respectively. We see that these distances are relatively 
small and decrease with increasing dimension for 2 < d < 5 and would expect the same 
trend to continue beyond five dimensions. This behavior is expected because RSA packings 
have superexponential decay of large-distance pair correlations in any dimension, and the 
decorrelation principle dictates that all unconstrained correlations at any pair distance must 
vanish as d becomes large. Therefore, a saturated RSA packing in high dimensions should 
be well approximated by an ensemble in which spheres randomly and sequentially packed in 
a local region around a centrally located sphere until saturation is achieved. Indeed, we have 
verified this proposition numerically, but do not present such results here. As the dimension 
increases, therefore, the local environment around a typical sphere in an actual RSA packing 
in IR d should be closely approximated by saturated RSA configurations on the surface of a 
(i-dimensional sphere, even though the former cannot have contacting particles. 



19 




FIG. 5: (Color online) The cumulative coordination number Z(r) for RSA packings for the first six 
space dimensions very near their respective saturation densities, as obtained from our numerical 
data for g2 (r) and (|14|) . The insets show the concavity of Z{r) near the contact value, which is in 
agreement with the behavior predicted by relation ()23|) . Away from contact, formula (|24jl provides 
a good approximation to the numerical data, especially for d > 3. 

D. Structure Factor 

In light of the discussion given in Section II, the total correlation function h(r) decays 
to zero for large r super-exponentially fast. It is well known from Fourier transform theory 
that if a real-space radial function f(r) in M d decreases sufficiently rapidly to zero for large 
r such that all of its even moments exist, then its Fourier transform f(k) is an even function 



20 



and analytic at k = 0. Thus, the structure factor S(k) for an RSA packing of spheres in R d 
must be an even, analytic function at k — 0. Hence, S(k), defined by (|16p. has an expansion 
about = in any space dimension d for < < <p s of the general form 

S(k) = S + S 2 k 2 + O(k 4 ), (25) 

where Sq and Si are the <i-dependent constants defined by 

POO 

So = 1 + 2 d d(j) / r d - 1 /i(r)rfr > (26) 



and 



oo 



S 2 = -2 d ~V \ r d+1 h{r)dr. (27) 



This analytic behavior of S(k) is to be contrasted with that of sphere packings near the 
MRJ state, which possesses a structure factor that is nonanalytic at k — 33| due to a 
total correlation function h(r) having a power-law tail. 

It is of interest to determine whether RSA packings are hyperuniform [23[ as <ft ~ * <Ps and, 
if not, their "distance" from hyperuniformity. A hyperuniform packing is one in which 

lim5(fc) -f 0, (28) 

i.e., the infinite- wavelength density fluctuations vanish. For RSA packings, this is equivalent 
to asking whether the generally nonnegative coefficient Sq, defined in (J25|) . vanishes. It is 
known that in one dimension, So ~ 0.05 and hence RSA rods are nearly but not quite 
hyperuniform. For any non-hyperuniform packing, the magnitude of So provides a measure 
of its "distance" from hyperuniformity. For a Poisson point pattern, it is well known that 
S = 1, but, in general, So can become unbounded if h{r) decays to zero more slowly than 
r~ d , as it does for a fluid at its critical point. 

Our results for S(k) very near the saturation density for the first six space dimensions 
are depicted in Fig. |3 To our knowledge, these results for d > 2 have not been presented 
before. As the space dimension increases, the amplitudes of the oscillations of S(k) diminish, 
consistent with the decorrelation principle. Note that the minimum value of S(k) in each 
dimension is achieved at the origin. Table ITTT1 provides the value of the structure factor at 
k = 0, denoted by So, by extrapolating our numerical data from the "direct" method near 
k = using the form ()25|) up to quadratic terms. We see that for 1 < d < 6, all of the 
packings are nearly hyperuniform and for d > 2, the distance from hyperuniformity does not 

21 



4 8 12 16 20 

kD 

2 
1.5 
^ 1 
0.5 

°0 4 8 12 16 20 

kD 

FIG. 6: (Color online) The structure factor S(k) for RSA packings for the first six space dimensions 
very near their respective saturation densities. Top panel includes curves for d = 1,2 and 3 and 
the bottom panel includes curves for d = 4, 5 and 6. Consistent with the behavior of we see 
again that pair correlations clearly decrease as the space dimension increases. 

appreciably vary as a function of dimension. In fact, our results indicate that the minimum 
value So quickly approaches a constant value of about 1/20 = 0.05 as d becomes large. 

The near hyperuniformity of a RSA packing in W 1 at its respective maximal density is 
a consequence of the saturation property. Long wavelength density fluctuations are appre- 
ciably suppressed because spherical gaps of a diameter equal to a sphere diameter or larger 
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TABLE III: The structure factor So at k = at the stopping density (j) s top obtained extrapolating 
our numerical data from the "direct" method near k = using the form (|25[) up to quadratic terms. 



Dimension, d 




1 


0.051 


2 


0.059 


3 


0.050 


4 


0.050 


5 


0.050 


6 


0.050 



cannot exist in the packing. In an earlier paper |23j, two of us (S. T. and F. H. S.) conjectured 
that saturated but strictly jammed disordered packings must be hyperuniform, which was 
subsequently verified by numerical simulations in three dimensions |23j. By this reasoning, 
one would expect a ghost RSA packing in R d not to be hyperuniform, even at its maximal 
density of l/2 d because it is never saturated. In fact, in Appendix B we determine the 
structure factor of the ghost RSA packing exactly and show that at its maximal density its 
distance from hyperuniformity increases as the space dimension d increases, asymptotically 
approaching the value of 1/2. 



V. HIGH-DIMENSIONAL SCALING OF SATURATION DENSITY AND LOWER 
BOUNDS FOR ANY d 

To determine whether the form of the observed density scaling (|2(Jj) for the first six space 
dimensions persists in the high-dimensional limit, we will apply an optimization procedure 



that we introduced to study the structure of disordered packings [8|, [14(. We begin by briefly 
reviewing this procedure and then apply it to the problem at hand. A g 2 -invariant process is 
one in which a given nonnegative pair correlation ^(r) function remains invariant as density 
varies for all r over the range of densities 

< < 0*. (29) 
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The terminal density 0* is the maximum achievable density for the (72-invariant process 
subject to satisfaction of certain necessary conditions on the pair correlation function. In 
particular, we considered those "test" g2(x) , s that are distributions on R d depending only 
on the radial distance x. For any test g2{x), we want to maximize the corresponding density 
4> satisfying the following three conditions: 

(i) 92(f) > for all r, 

(ii) 92(f) = for r < D, 
(iii) 

d f°° JdAkr) 
S(k) = l + p (2?r) 5 / r d - x h(r) 2 A _ dr > for all k, 
Jo (At) 2 1 

where S(k) is the structure factor defined by (|15j) . When there exist sphere packings with 
92 satisfying conditions (i)-(iii) for </> in the interval [0, 0*], then we have the lower bound on 
the maximal density given by 

0max > 4>*- (30) 

In addition, to the nonnegativity of the structure factor S(k), there are generally many 
other conditions that a pair correlation function of a point process must obey [35] . There- 
fore, any test 92 that satisfies the conditions (i)-(iii) does not necessarily correspond to a 
packing. However, two of us (S. T. and F. H. S.) have conjectured that a hard-core nonneg- 
ative tempered distribution q 2 (r) is a pair correlation function of a translationally invariant 
disordered sphere packing |l9j ] in M. d at number density p for sufficiently large d if and only 
if S(k) > The maximum achievable density is the terminal density 0*. 

A certain test 92 and this conjecture led to the putative long-sought exponential improve- 
ment on Minkowski's lower bound The validity of this conjecture is supported by a 
number of results. First, the decorrelation principle states that unconstrained correlations 
in disordered sphere packings vanish asymptotically in high dimensions and that the g n for 
any n > 3 can be inferred entirely from a knowledge of p and g 2 - Second, the necessary 
Yamada condition [36] appears to only have relevance in very low dimensions. This states 
that the variance a 2 (Q) = ((N(Q) 2 — (N(Q))) 2 ) in the number N(Q) of particle centers 
contained within a region or "window" Q C M. d must obey the following condition: 

a 2 (tt) = p\tt\ U + P [ h(r)dr] > 9(1 - 9), (31) 

where 9 is the fractional part of the expected number of points p\Q\ contained in the window. 
Third, we have shown that other new necessary conditions also seem to be germane only 
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in very low dimensions. Fourth, we have recovered the form of known rigorous bounds on 
the density in special cases of the test #2 when the aforementioned conjecture is invoked. 
Finally, in these latter two instances, configurations of disordered sphere packings on the 
flat torus have been numerically constructed with such g 2 in low dimensions for densities up 



to the terminal density 
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We now apply this optimization procedure and the aforementioned conjecture to ascertain 
whether the form of the density scaling (j20j) persists in the high-dimensional limit. In this 
limit, the decorrelation principle as well as our results for the pair correlation function in the 
first six space dimensions enable us to conclude that g2{x) is very nearly unity for almost all 
distances beyond contact except for a very small nonnegative interval in the near-contact 
region. Therefore, because the extra structure in low dimensions representing unconstrained 
spatial correlations beyond a single sphere diameter should vanish as d —>■ oo, we consider a 
high-dimensional test pair correlation function in IR d that is nonunity within a small positive 
interval 1 < x < 1 + e beyond contact and unity for all x greater than 1 + e, i.e., we consider 



92 [x) 



0, x < 1, 

l + f(x), l<x<l + e, (32) 

1, x > 1, 



where e is a nonnegative constant (e > 0) and f(x) is any integrable function in one dimen- 
sion that satisfies f(x) > — 1. This class of functions can include even those that diverge to 
infinity as x — > 1. Examples of the latter integrable class include 

/(x) = -lm>-l), (33) 

fix) = < a < 1, (34) 

and 

f(x) = S(x-l), (35) 

where S(x) is the Dirac delta function. Equation ()33|) describes the divergence seen in gi{x) 
of the standard RSA packing at contact [cf. (JSj)]. By contrast, Eq. (J53jl characterizes the 
near-contact divergence of g2{x) for maximally random jammed (MRJ) packings [sj] with 

o.6 QQ- 

Equation (jHHj) describes random sphere packings with a positive average 
contact number. For general f(x), the corresponding structure factor [cf. (iii)] for the test 
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function (|32|) in any dimension d is given by 

2¥0r(i + d/2) 



S(k) = 1 



k 2 



1+e 



Jd(k) — k I x 2 f(x)Jd_ l (kx)dx 



(36) 



where i/ = d/2. 

We now make use of the following general result that applies to any function G(x) that 
is bounded clS X ^ X ! 



lim 



1+e 



where I(x) is the indefinite integral 



G(x)f(x)dx = G(l)I(e) 



I(x) = / f(x)dx. 



For the functions (|3^j) . (|3*4*|) and (|33j). the integral /(x = e) is respectively given by 



J(e) = (l-lne)e, 



1 — a ' 



< a < 1. 



and 



J(x) = 1. 

Making use of the result (pT7jl in (jHEJ) yields the structure factor to be given by 



1 2¥0r(i + d/2) 



Jd(k) - Ud^ik)!^) 



The structure factor for small k can be expanded in a MacLaurin series as follows: 

S(k) = 1 + 2 d <j)[dl(e) - 1] + - I(e)(d+ 2)] /c 2 + 0(fc 4 ). 

The last term changes sign if 1(e) increases past l/(d + 2). At this crossover point, 



S(k) = 1 - 



2 d+i 
d + 2 



+ 0(k 4 



(37) 



(38) 



(39) 
(40) 

(41) 
(42) 



(43) 



(44) 



Under the constraint that the minimum of S(k) occurs at k = 0, the terminal density is 
then given by 



d + 2 



:i - So) 



(45) 



where So G [0, 1] is the value of the structure factor at k = or the assumed minimum 
value in the high- dimensional limit. Thus, we see that the terminal density is independent 
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of the specific form of /(e) or, equivalently, the choice of the function /(r) [cf. (J3~2*|) ]. which 
only has influence in a very small nonnegative interval around contact. For a hyperuniform 
situation (So = 0), the formula (f£f|) reduces to 

*. = ^ w 

which was obtained previously for the specific choice of /(r) given by (|35|) . 

It is noteworthy that the high- dimensional asymptotic structure factor relation (|42jl under 
the conditions leading to (j^S) yields a structure factor for d = 6, a relatively low dimen- 
sion, that is remarkably close to our corresponding simulational RSA result (depicted in the 
bottom panel of Fig. [HJ for most values of the wavenumber k. Such agreement between 
the asymptotic and low-dimensional numerical results strongly suggests that our asymp- 
totic form (|32j) for the pair correlation function indeed captures the true high-dimensional 
behavior for RSA packings. 

In summary, we see that the high-dimensional result (f4"5|) shows that the form of the 
density scaling (J2Uj) at saturation for relatively low- dimensional RSA packings persists in 
the high- dimensional limit. Indeed, as we will show below, the high-dimensional scaling ()45|) 
provides a lower bound on the RSA saturation density <p s for arbitrary d, i.e., 

*.>^(l-So). (47) 

Figure [7| depicts a graphical comparison of the lower bound (|4"7jl to our numerical data for 
the saturation density for the first six space dimensions. 

The lower bound (|4"7|) is a consequence of a more general principle that enables us to 
exploit high-dimensional information in order to infer scaling behavior in low dimensions, as 
we now describe. For any particular packing construction in M. d (e.g., RSA, ghost RSA or 
MRJ packings), the highest achievable density <pm{d) decreases with increasing dimension 
d. Therefore, the scaling for the maximal density in the asymptotic limit d — > oo, which we 
denote by m (oo), provides a lower bound on 4> m (d) for any finite dimension, i.e., 

(j) m (d) > m (oo). (48) 

In the case of RSA packings in M d , we have shown that the high-dimensional density 
scaling is provided by the analysis leading to flj5)) and therefore use of (|4*8Jl yields the lower 
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FIG. 7: (Color online) Comparison of the lower bound (|47|1 on the saturation density (with Sq = 
0.05) to our corresponding numerical data for the first six space dimensions. 

bound (}4"Tj) in which So is a small positive number. It is noteworthy that a lower bound on 
the MRJ density 4>mrj is given by the right side of inequality of (}4*Tj) but with So = 0, i.e., 

d + 2 



'MRJ 



> 



2 d+l 



(49) 



This is obtained by recognizing that the same high- dimensional scaling analysis as we used 
for RSA applies with one qualitative difference. The high- dimensional limit of the pair 
correlation function of an MRJ packing is expected to be of the same form as (}32|) but where 
fix) is a Dirac delta function to account for interparticle contacts due to the constraint of 

nn 

jamming. Since we know that MRJ packings are hyperuniform (So = 0) |33j,|3J], then (|45|) 
together with (gHJ) produces the bound flUj). The MRJ lower bound (JH yields 0.3125, 
0.1875, 0.109375, 0.0625 for d = 3,4,5 and 6, respectively, which is to be compared to the 



corresponding actual MRJ densities of 0.64, 0.46, 0.31 and 0.20 |31|, |32j, |33j, |3J . Note that 
in Ref. JjJ, the right side of ()49|) was argued to be a lower bound on the maximal density 
</Wx of any sphere packing in M. d . In the case, of the ghost RSA packing, we know that l/2 d 
is the maximal density <Pgrsa for any dimension, and therefore this result in conjunction 
with (J4*%j) yields the lower bound <J)grsa > l/2 d , which of course is exact. 
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VI. CONCLUSIONS 



We have studied the structural characteristics of random sequential addition (RSA) of 
congruent spheres in (^-dimensional Euclidean space M. d in the infinite-time or saturation limit 
for the first six space dimensions (1 < d < 6) both numerically and theoretically. Specifi- 
cally, by numerically generating saturated RSA configurations in each of these dimensions, 
we determined the saturation density, pair correlation function, cumulative coordination 
number and the structure factor. We found that for 2 < d < 6, the saturation density 
(f) s has the scaling given by (f2*U|) . Using theoretical considerations, we showed analytically 
that the same density- scaling form persists in the high-dimensional limit. A byproduct of 
the aforementioned high- dimensional analysis was the determination of a relatively sharp 
lower bound on the saturation density ()47|) of RSA packings for any d, which utilized the 
infinite- wavelength limit of the structure factor in the high- dimensional limit. Thus, high- 
dimensional information was exploited to provide density estimates in low dimensions. We 
proved rigorously that a Palasti-like conjecture cannot be true for RSA hyperspheres. We 
also demonstrated that the structure factor S(k) must be analytic at k = and that RSA 
packings for 1 < d < 6 are nearly "hyperuniform" (i.e., infinite wavelength density fluc- 
tuations vanish). Consistent with the recent "decorrelation principle," we find that pair 
correlations markedly diminish as the space dimension increases up to six. 

In Appendix A, we obtained kissing number statistics for saturated RSA configurations on 
the surface of a <i-dimensional sphere for dimensions 2 < d < 5 and compared to the maximal 
kissing numbers in these dimensions. The discrepancy between average RSA kissing numbers 
and maximal kissing number was found to increase as the space dimension increased. Finally, 
in Appendix B, we determined the structure factor exactly for the related "ghost" RSA 
packing in M. d and showed that its distance from "hyper uniformity" increases as the space 
dimension increases, approaching a constant asymptotic value of 1/2. 

It is interesting to observe that the best known rigorous lower bound on the maximal 
density [3], derived by considering Bravais lattice packings, has the same form as the den- 
sity scaling (|2U|) for RSA packings, i.e., for large d, it is dominated by the term d/2 d . The 
fact that the saturation density of disordered RSA packings approaches this rigorous lower 
bound suggests the existence disordered packings whose density surpasses the densest lat- 
tice packings in some sufficiently high dimension. The reason for this is that we know that 



29 



there are disordered packings in low dimensions whose density exceeds that of correspond- 
ing saturated RSA packings in these dimensions, such as maximally random jammed (MRJ) 



packings 
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33 



34j . The density of saturated RSA packings in dimension d is sub- 



stantially smaller than the corresponding MRJ value because, unlike the latter packing, the 
particles can neither rearrange nor jam. The possibility that disordered packings in suffi- 
ciently high dimensions are the densest is consistent with a recent conjectural lower bound 
on the density of disordered hard-sphere packings that was employed to provide the putative 
exponential improvement on Minkowski's 100-year-old bound 14]. The asymptotic behavior 
of the conjectural lower bound is controlled by 2~(°- 77865 ---)' i . 

Challenging problems worth pursuing in future work are the determinations of analytical 
constructions of disordered sphere packings with densities that equal or exceed d/2 d for 
sufficiently large d or, better yet, provide exponential improvement on Minkowski's lower 
bound. The latter possibility would add to the growing evidence that disordered packings 
at and beyond some sufficiently large critical dimension might be the densest among all 
packings. This scenario would imply the counterintuitive existence of disordered classical 
ground states for some continuous potentials in such dimensions. 
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APPENDIX A: RSA KISSING NUMBER 

The number of unit spheres in M. d that simultaneously touch another unit sphere without 
overlap is called the kissing number (also known as the contact number or coordination 
number). The kissing number problem seeks the maximal kissing number Z max as a function 
of d. In dimensions 1 < d < 3, the maximal kissing numbers are known and correspond to 
the kissmg nu m b erS ofVdense, sphere pa^, wh.ch a r e Bravais lattices fl Q: the 
simple linear lattice for d — 1, the triangular lattice for d = 2, and the face-centered cubic 
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(FCC) lattice for d = 3. For d = 1 and d = 2, the maximal kissing numbers are exactly 2 and 
6, respectively. Although one dimension is trivial, it is a "miracle" of two dimensions that 
6 circles can simultaneously touch another circle without any gaps, and in this sense this 
unique configuration (up to trivial rotations) is "rigid" because there are no displacements of 
the 6 contacting circles that lead to a different configuration while maintaining the contacts. 
This unique kissing number arrangement is also six-fold symmetric. In three dimensions, it 
is known that Z max = 12, which is achieved by the FCC sphere packing, but there are no 
unique configurations because gaps exist between contacting spheres that enable one optimal 
kissing configuration to be displaced into a different optimal configuration, and therefore 
optimal configurations need not have any symmetry. The aforementioned subtleties in the 
three-dimensional case was at the heart of a famous debate in 1694 between Issac Newton 
(who claimed that Z max = 12) and David Gregory (who contended that Z max = 13). 

One of the generalizations of the FCC lattice to higher dimensions is the checkerboard 
lattice, defined by taking a cubic lattice and placing spheres on every site at which the sum of 
the lattice indices is even i.e., every other site. The densest packing for d — 4 is conjectured 
to be the -D4 lattice, with a kissing number Z = Z max = 24 [5], which is also the maximal 
kissing number in d — 4 |4l| . This optimal configuration is referred to as the 24-cell, which is 
both rigid and highly symmetric. For d = 5, the densest packing is conjectured to be the D§ 
lattice with kissing number Z = 40 [5|. This kissing configuration is also highly symmetrical. 
The maximal kissing numbers Z max for d — 5 is not known, but has the following bounds: 
40 < Z max < 46. 

Here we determine the distribution of kissing numbers Z; by placing hyperspheres ran- 
domly and sequentially on the surface of a hypersphere at the origin until the surface is 
saturated. We call such a configuration a saturated RSA kissing number configuration. The 
average kissing number (Z) is given by 

(Z) = J2z t P(Z t ), (Al) 
1=1 

where P{Zi) the probability of finding a saturation kissing number Zj. We begin our sim- 
ulations by placing a central hypersphere of unit diameter at the origin of a hypercubic 
simulation box. A large number of points n pts (n pts = 2 x 10 5 for d = 2, n pts = 10 6 for 
2 < d < 4, and 10 7 points for d = 5) are uniformly distributed in the region between the 
exclusion hypersphere of unit radius surrounding the hypersphere of unit diameter and the 
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boundary of the simulation box. Each point is randomly and sequentially radially projected 
(in the direction of the hypersphere center) to the surface of the exclusion hypersphere (via 
a radial distance rescaling) subject to the nonoverlap condition, i.e., a projected point is 
accepted if the angular separation between it and any other previously accepted point is 
greater than or equal to 60 degrees, otherwise it is rejected. The simulation terminates 
when all projected points obey this nonoverlap condition at which time the exclusion-sphere 
surface is taken to be saturated, i.e., the surface of the central sphere of unit diameter is 
saturated with contacting spheres of unit diameter. We found that the number of points n pts 
that we initially distributed in the simulation box before the projection step is sufficiently 
large to ensure that the surface of the hypersphere is truly saturated after the projection 
step in the dimensions considered. 

Table ITVl provides kissing number statistics for saturated RSA configurations for dimen- 
sions 2 < d < 5. For d — 2, only two values of Zj are allowed: 4 and 5. Any kissing 
number equal to 3 or less is prohibited because such configurations are not saturated. On 
the other hand, a kissing number of 6 has a probability of zero of occurring by random 
sequential addition because the configuration corresponding to this optimal kissing num- 
ber is unique. We find that the average kissing number is approximately equal to 4.5. For 
d = 3, the average kissing number (Z) = 8.34135. The fact that the maximal kissing number 
configurations (Z max = 12) in M. 3 are nonunique implies that a configuration of 12 spheres 
has a positive (albeit small) probability of occurring via an RSA process. Nonetheless, we 
were not able to observe such a configuration in a total of 10 6 configurations. The smallest 
observed kissing number for d = 3 was six, which presumably is the smallest number re- 
quired for a saturated packing. The minimal kissing number configurations for saturation 
are related to the following problem: How can n points be distributed on a unit sphere such 
that they maximize the minimum distance between any pair of points? For six points, the 
solution to his problem is well known: they should be placed at the vertices of an inscribed 
regular octahedron. Since the minimum angular separation between any pair of points in 
this highly symmetric case is 90 degrees, the associated kissing number configuration is sat- 
urated. Note that highest kissing number of 18 reported for d = 4 is substantially smaller 
than the maximal kissing number Z mSuX = 24. Apparently, achieving kissing numbers that 
approach those of the optimal highly symmetric, rigid 24-cell configuration by a random 
sequential addition is effectively impossible. Therefore, it is plausible that the possible con- 
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TABLE IV: Kissing number statistics for saturated RSA configurations on the surface of a d- 
dimensional sphere for dimensions 2 < d < 5. Here Zj is the integer- valued saturation kissing 
number, -P(Zj) the probability of finding a saturation kissing number Zj, (Z) the average saturation 
kissing number. In each dimension, the statistics are determined from 100,000 configurations. We 
also include the largest known kissing numbers Z max . 



Dimension, d 


Kissing Number, Zi 


Probability, P(Z;) 


(Z) 


Zmax 


2 


4 


0.515680 


4.48432 


6 




5 


0.484320 






3 


6 


0.001400 


8.34957 


12 




7 


0.091020 








8 


0.502230 








9 


0.367400 








10 


0.037860 








11 


0.000090 






4 


11 


0.001840 


13.80530 


24 




12 


0.055960 








13 


0.302440 








14 


0.435910 








15 


0.183060 








16 


0.020260 








17 


0.000520 








18 


0.000010 






5 


17 


0.000030 


21.46765 


40 




18 


0.001530 








19 


0.024510 








20 


0.146930 








21 


0.341250 








22 


0.329250 








23 


0.132350 








24 


0.022290 








25 


0.001810 








26 


0.000050 







figurations corresponding to kissing numbers of 19-23 are also characterized by high degree 
of symmetry (and possibly rigidity) based upon the absence of such kissing numbers. The 
fact that the average kissing number for d = 5 is substantially lower than the highest known 
kissing number of 40 is presumably related to the high symmetries required to achieve high 
Z values in this dimension. 

Our data for the average RSA kissing number over the range of considered dimensions is 
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d 

FIG. 8: (Color online) Our numerical data (black circles) for the average kissing number (Z) as a 
function of dimension d and the quadratic fit function ()A2|) (solid curve). 

fit very well by the following quadratic expression in d: 

(Z) = b + b 1 d + b 2 d 2 } (A2) 

where 60 = 2.74488, b\ = —1.01354 and b 2 = 0.950395, and the correlation coefficient is 
0.9999. The data and this quadratic fit function are depicted in Figure |H1 If this expression 
persisted for large d, it would predict that the average RSA kissing number asymptotically 
grows as d 2 . 



APPENDIX B: STRUCTURE FACTOR FOR GHOST RSA PACKINGS 

For ghost RSA packings of spheres of diameter D in W 1 , the pair correlation function in 
;he infinite-time limit is given exactly for any space dimension d by the following expression 

Q(r-D) 
1 - a 2 {r;D)/2 

where Q(x) is the unit step function, equal to zero for x < and unity for x > 1, and 
a 2 (r; D) is the intersection volume of two spheres of radius D whose centers are separated 
by the distance r divided by the volume of a sphere of radius D. Expressions for the scaled 
intersection volume a 2 (r; D) for any d are known exactly; see Refs. and for two 
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different representations. The scaled intersection volume a 2 {r; D) takes its maximum value 
of unity at r = and monotonically decreases with increasing r such that it is nonzero 
for < r < 2D, i.e., it has compact support. The corresponding total correlation h(r) = 
g 2 {r) — 1 is given by 

h(r) = -Q(D - r) + OL 2 (r; D) i _ D ^ 2D _ ^ (B2) 

2 — a 2 (r; D) 

We see that that h(r) can be written as a sum of two contributions: the pure step function 
contribution — Q(D — r), which has support for < r < D, and a contribution involving 
a 2 (r; D), which has support D < r < 2D. Substitution of (jB2|) into (fTBj) yields the structure 
factor to be 

S(k) = S SF (k) + S EX (k), (B3) 

where 

Jd(kD) 

S SF (k) = l-2^T(l + d/2)^— - (B4) 

(kD)2 

is the structure factor for the step function contribution —Q(D — r), and 

S EX (k) = 2iT( l + m fr'-i ">< r f Ji '-f\ r (B5) 

Jd 2-a 2 {r;D) (kr)?' 1 

is the contribution to S(k) in excess to the structure factor for the step function. Here 
we have used the fact that the infinite-time density is <p — l/2 d . For odd dimensions, 
SEx{k) can be obtained explicitly in terms of sine, cosine, sine integral and cosine integral 
functions. We do not explicitly present these expressions here but instead plot S(k), defined 
by (|B3|) . for various dimensions in Fig. For d — 1,3 and 11, S(k = 0) is given by 
0.150728, 0.290134 and 0.452217, respectively, and therefore not only is the ghost RSA 
packing not hyperuniform, as expected, but its distance from hyperuniformity increases 
as the space dimension d increases, asymptotically approaching the value of 1/2. This 
should be contrasted with the standard RSA packing in M. d , which we have shown is nearly 
hyperuniform. 

In the limit d — > oo, the excess contribution to the structure factor has the limiting form 



S E x(k) -> - 



2?r(l + d/2)Ji(kD) 



(kD) 



2 



S SF (k)] 2 . (B6) 
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FIG. 9: (Color online) The structure factor S(k) versus kD for various space dimensions (d = 1,3 
and 11) for ghost RSA packings. 

The resulting structure factor in this asymptotic limit is given by 

S(k) = S SF (k) + S EX (k) 

= \[l + S 2 SF (k)}, d^oo. (B7) 

It has the following small-fc expansion: 

s ^-\ + WTW k '- nd + 2n d+i f + ^ k ^°- < B8) 

The asymptotic result ()B6|) is easily obtained by utilizing the fact that in the limit d — > oo, 
a 2 (r; _D)/[2 — a 2 (r; D)] — > a 2 (r; D)/2 jl^j]. Substitution of this result into the general relation 
()B5j) and recognizing that the lower limit D of this integral can be replaced by in the limit 
d — > oo yields the asymptotic form 

a 2 (k;D) 

SExik) "^pp (B9) 

where a 2 (fc; D) denotes the Fourier transform of a 2 (r; D) and V\{D) is the volume of a sphere 
of radius D [cf. (J2J)]. The quantity a 2 (k;D) is known explicitly in any dimension [23 1 and 
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substitution of this result into (jB9|) immediately yields (jB6|) . 
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